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ABSTRACT 

We present new FUV observations of the pulsating DA white dwarf G226— 29 obtained with the 
Far Ultraviolet Spectroscopic Explorer (FUSE). This ZZ Ceti star is the brightest one of its class and 
the coolest white dwarf observed by FUSE. We report the first detection of the broad quasi-molecular 
collision-induced satellite of Ly-/3 at 1150 A, an absorption feature that is due to transitions which take 
place during close collisions of hydrogen atoms. The physical interpretation of this feature is based on 
recent progress of the line broadening theory of the far wing of Ly-/3. This predicted feature had never 
been observed before, even in laboratory spectra. 

Subject headings: line: profiles - stars: individual (G226— 29) - stars: white dwarfs - stars: 
atmospheres - ultraviolet: stars 



1. INTRODUCTION 

DA white dwarfs have hydrogen-rich atmospheres whose 
FUV spectra show the Lyman series lines. Spectra of 
these high gravity stars are modeled by theoretical spectra 
based on the Stark broadening of the hydrogen lines (see, 
e.g., Vidal, Cooper, & Smith 1973; Finley et al. 1997; 
Barstow et al. 2003). FUV and UV observations of some 
DA white dwarfs show, however, strong deviations from 
the Stark broadening. Greenstein (1980) and Holm et 
al. (1985) reported that International Ultraviolet Explorer 
(IUE) spectra of cool white dwarfs show strong absorption 
features in the far wing of Ly-a near 1400 A and 1600 A. 
Koester et al. (1996) announced the presence of similar 
features in the red wing of Ly-/3 by analyzing observations 
of Wolf 1346 obtained with the Hopkins Ultraviolet Tele- 
scope. More recently, absorption features in the wing of 
Ly-7 were detected in the Orbiting and Retrievable Far 
and Extreme Ultraviolet Spectrometers (ORFEUS) and 
Far Ultraviolet Spectroscopic Explorer (FUSE) spectra of 
white dwarfs (Koester et al. 1998; Wolff et al. 2001; 
Hebrard et al. 2003; Dupuis et al. 2003). All these absorp- 
tion features are interpreted as quasi-molecular satellites 
of Ly-a (Koester et al. 1985; Nelan & Wegner 1985), Ly- 
P (Koester et al. 1996), or Ly-7 (Hebrard et al. 2003). 
Quasi-molecular lines arise from radiative collisions of ex- 
cited atomic hydrogen with unexcited neutral hydrogen 
atoms or protons. 

Theoretical calculations of the complete Ly-/3 line pro- 
file which include perturbations by both neutral hydrogen 
and protons (Allard et al. 2000) have been used recently to 
improve theoretical modeling of synthetic spectra for cool 
DA white dwarfs (Hebrard et al. 2002a). These new cal- 
culations reveal in particular a broad H2 collision-induced 



satellite in the red wing of Ly-/3 at 1150 A. This satellite 
appears in models covering a narrow range of effective tem- 
peratures and is strong enough to be detected in spectra of 
cool DA white dwarfs with 10000 < T cS < 13000 K. In ad- 
dition to that, the H2 A1150 satellite is very sensitive to the 
degree of ionization and represents a potentially important 
means to diagnose the convective mixing efficiency in DA 
white dwarfs, especially in this range of effective tempera- 
tures (see, e.g., Bergeron et al. 1992). Even in laboratory 
work, where a laser-produced plasma generates controlled 
conditions similar to those encountered in stellar atmo- 
spheres of cool white dwarfs, this transition has yet to be 
detected. 

The effective temperature domain in which the H2 A1150 
satellite is the strongest includes an important class of 
stars: the ZZ Ceti stars. ZZ Ceti stars are a class of 
hydrogen-rich atmosphere white dwarfs exhibiting multi- 
periodic variations with period ranging from 100 to 1000 s 
that are produced by non-radial g-modes (see, e.g., the 
recent review by Fontaine, Brassard, & Charpinet 2003). 
According to Bergeron et al. (1995) all ZZ Ceti stars oc- 
cupy an empirical instability strip delimited by 11100K < 
T e ff < 12500K. These objects have a great astrophysical 
interest because seismological studies allow to probe their 
interiors and therefore provide insights for constraining the 
stellar structure and evolution of white dwarfs. ZZ Ceti 
stars represent an evolutionary phase through which all 
DA stars must evolve. Hebrard et al. (2002a) have pre- 
dicted that both H2 and Hj satellites of the Ly-a and 
Ly-/3 lines should be detectable in ZZ Ceti stars. 

Far-UV spectra have never been explored for any 
ZZ Ceti star. The observations obtained with IUE and 
HST were limited to the red wing of Ly-a. Here, we 
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present a far-UV spectrum (down to 1000 A) of the bright- 
est object of this class, G226— 29. These data allow us to 
report the first detection of the H 2 Ly-/3 collision-induced 
(CI) satellite at 1150 A. 

2. THE COLLISION-INDUCED SATELLITE AT 1150 A 

Based on a theoretical study of the Ly-/3 profile of 
atomic hydrogen perturbed by collisions with neutral hy- 
drogen atoms and protons (Allard et al. 2000), we pre- 
dict a CI satellite feature of H2 at 1150 A. In summary, 
our theoretical approach is based on the quantum theory 
of spectral line shapes of Baranger (1958a, 1958b) devel- 
oped in an adiabatic representation to include the degen- 
eracy of atomic levels (Royer 1974; Royer 1980; Allard 
et al. 1994). Key ingredients in our calculations are the 
potential energies V(R) for each electronic state of the 
or H 2 molecule (R denotes the internuclear distance 
between the radiator and the perturber). For H-H + colli- 
sions, we have used the potentials calculated by Mad- 
sen & Peek (1971) and, for H-H, the H 2 potentials cal- 
culated by Detmer et al. (1998) and Schmelcher (2000). 
The resulting complete Ly-/3 profile is shown in Fig. 1. 
The broad line satellite near 1150 A is due to the pertur- 
bation in an atomic hydrogen collision corresponding to 
the B"B *E+ - X X S+ transition of H 2 . 

The properties of the H 2 A1150 CI satellite can be bet- 
ter understood by studying Fig. 2, which illustrates the 
difference potential energies AV(R) for this transition: 

AV(R) = V e , e (R) = V e ,(R)-V e (R) , (1) 

where e and e' label the energy surfaces on which the inter- 
acting atoms approach the initial and final atomic states 
of the transition as R — > 00. The most significant charac- 
teristic of the potential curve is the existence of the double 
wells. Each extremum in the difference potential leads, in 
principle, to a corresponding satellite feature in the wing 
of Ly-/3. Another important feature is that at larger inter- 
nuclear separation, up to R = 19 A, the state has an ionic 
character. The ionic interaction weakens slowly (1/i? 2 ver- 
sus 1/R 6 dependence) thus making the potential energy 
difference broad compared to the steep well of the B-X 
transition which gives rise to the 1600 A satellite in the 
Ly-a wing. This different shape is important because the 
position of the extremum and the functional dependence 
of the potential difference on internuclear separation de- 
termine the amplitude and shape of the satellites (Allard 
et al. 1994). 

The satellite amplitude depends also on the value of the 
electric dipolc transition moment D(R) (dashed line on 
Fig. 2) taken between the initial and final states of the 
radiative transition. The dependence of the moment on 
the separation of the atoms during a collision modifies rel- 
ative contributions to the profile along the collision trajec- 
tory. Line profile calculations have been done in a classical 
path theory which takes into account the variation of the 
electric dipolc moment during a collision (sec, Allard et 
al. 1999). For the transitions contributing to Ly-a and 
Ly-/3, we used the dipole moments for and H 2 that 
Ramaker & Peek (1972) and Spielfiedel (2003) calculated 
as a function of internuclear distances. 

The qualitative effects of the B"B — X transition can 
be discerned in Fig. 2, which shows the radiative dipolc 



transition moment and the difference potential energies as 
a function of R. The transition dipole moment is small for 
most values of R but shows a broad maximum at 4.5 A, 
close to the broad minimum of the outer well in the po- 
tential difference. In contrast the dipole moment is small 
in the vicinity of the inner well in the potential differ- 
ence. As a result, the main contribution from this tran- 
sition is a CI satellite in the far wing centered at a fre- 
quency corresponding to the broad minimum (Allard and 
Kielkopf 1982). 

The CI absorption depends strongly on the internuclear 
separation and produces a broad spectral feature with a 
characteristic width of the order of the inverse of the du- 
ration of the close collision. 

The line satellite shown in Fig. 1 presents a shoulder at 
1120 A, a similar shape has been obtained for the 1600 A 
satellite. In Fig. 6 of Allard et al. (1999) both theory 
and experiment show an oscillatory structure between the 
satellite and the line, with a minimum at about 1525 A. 
These oscillations are an interference effect (Royer 1971; 
Sando & Wormhoudt 1973), and are expected to depend 
on the relative velocity of the collision and therefore on 
temperature. To conclude, we emphasize the importance 
of the accuracy of both the potential energies and the 
dipole moments for the line shape calculations and of a 
theory that takes into account the variation of the dipolc 
moment during an atomic collision. 

3. FUSE OBSERVATION OF G226-29 

G226-29 (WD1647+591, DN Dra) is the brightest 
ZZ Ccti star. It was originally discovered to be variable by 
McGraw & Fontaine (1980) but the first detailed analysis 
of the light curve had been published by Kepler, Robin- 
son, & Nathcr (1983). This is one of the best studied 
ZZ Ceti stars, for which Fontaine et al. (1992) proposed 
that pulsation properties could be interpreted assuming a 
"thick" hydrogen layer mass of logg(H) w —4.4. Subse- 
quent observations by the Whole Earth Telescope consor- 
tium (Kepler et al. 1995) and by HST/FOS (Kepler et al. 
2000) allowed unequivocal mode identification and sup- 
ported Fontaine et al.'s (1992) contention that G226— 29 
has a thick hydrogen layer. G226— 29 defines the blue (hot) 
edge of the instability strip according to optical determi- 
nation of Bergeron et al. (1995). 

G226-29 was observed by FUSE as part of our Cycle 4 
Guest Investigator Program D101. This ZZ Ceti star is the 
coolest white dwarf observed with this instrument so far. 
Two observations of four exposures each were obtained on 
2003 March 11 and 12 in time-tagged photon address mode 
(TTAG) with the object in the 30" x30" aperture (LWRS). 
The total duration was 29.1 ksec (~ 8 h). Details of the 
FUSE instrument may be found in Moos et al. (2000) and 
Sahnow et al. (2000).' 

The one-dimensional spectra were extracted from the 
two-dimensional detector images and calibrated using ver- 
sion 2.4.1 of the CalFUSE pipeline. The FUSE detec- 
tor segments of the two observations were co-added and 
projected on a 0.6 A-pixel base, i.e. pixels about 100 
times larger than the original FUSE detectors pixels. This 
degradation of the FUSE spectral resolution (typically 
A/AA ~ 15000 for this kind of target in the large slit; 
see Hebrard et al. 2002b; Wood et al. 2002) is of no effect 
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on the shapes of the large stellar features that we study 
and allows us to increase the signal-to-noise ratio. At this 
resolution, no spectral shifts were detected between the 
different co-added spectra. 

The two shorter wavelength segments (SiClB and 
SiC2A, ~ 905 A— 1000 A) were not used as no significant 
stellar spectrum was detected on them. We also did not 
use the LiFIB segment in our final spectrum, which is 
known to present a large-scale distortion (the worm) in the 
flux calibration that spans the region of interest (see, e.g., 
Hebrard et al. 2002a). Thus, the range 1105 A-1180 A, in 
which the H 2 satellite is detected (Sect. 4), comes from the 
LiF2A segment only. We checked that the H2 satellite is 
also detected on LiFIB, despite the alteration of the worm. 
This redundancy in the FUSE spectral coverage allow us 
to conclude that the feature around 1150 A is real, and 
not the result of an instrumental artifact. 

The final reduced spectrum is plotted in Fig. 3. 

4. SYNTHETIC SPECTRA AND COMPARISON WITH THE 
OBSERVATION 

We computed a small grid of LTE model atmospheres 
with pure hydrogen composition that explicitly include the 
Ly-a and Ly-/3 quasi-molecular opacities. These are deter- 
mined for T = 12000 K, computed from absorption profiles 
that take into account a variable dipole, modulated by the 
Boltzmann factor (Allard et al. 1999). The atmosphere 
models have been calculated using the program TLUSTY 
(Hubeny 1988; Hubeny & Lanz 1992, 1995). The syn- 
thetic spectra were computed by using the spectral syn- 
thesis code SYNSPEC. The line satellites have a strong 
blanketing effect and have been included in both the atmo- 
sphere model and synthetic spectra calculations. Convec- 
tion was treated within the usual framework of the mixing 
length approximation (MLT), where we used ML2/a=0.6 
(Bergeron et al. 1995). 

Fig. 3 compares our best model to the FUSE spectrum 
of G226-29. The observed feature at 1150 A is well re- 
produced by the predicted H 2 satellite, both in terms of 
shape and wavelength position. This is the first detection 
and identification of that feature. We also computed a 
model without including the H-H interactions to appre- 
ciate how these interactions affect the far wing of Ly-/3. 
Fig. 3 shows that this model (dotted line) neither pro- 
duces the H 2 A1150 absorption feature nor the significant 
flux absorption down to 1085 A. This illustrates that an ac- 
curate fit of the observed FUV spectra cannot be obtained 
without using the improved theoretical profile calculations 
including perturbations by both protons and neutral hy- 
drogen presented by Allard et al. (2000). 

As expected in that temperature range, the H^ satellite 
at 1076 A is also detected in spite of the low flux level. 
The satellite at 1058 A is not clearly detected because 
of the too low signal-to-noise ratio. The H~ absorption at 
1130 A (Wishart 1979) is detected; that H~ feature was 
already detected in the FUSE spectrum of the cool DA 
white dwarf G 231-40 (Hebrard et al. 2002a). 

The best fit is obtained for T cS = 12000 K and log g = 
7.9. That temperature is slightly lower than the one that 
Bergeron et al. (1995) determined from optical obser- 
vations (T e ff = 12460 ± 230 K), and it is similar to the 
one reported by Koester & Holberg (2001), who fitted the 



HST/FOS data by adding the constraints of the V mag- 
nitude and the trigonometric parallax (T c s— 12020 K). 
Our gravity, however, is significantly lower than those ob- 
tained in both studies (log g = 8.2 — 8.3). The causes of 
this difference arc not fully understood yet. We are now 
investigating these causes and hope to resolve the discrep- 
ancy between our value of the gravity and those of Berg- 
eron et al. (1995) and Koester & Holberg (2001). This 
work is beyond the scope of the present paper and will 
be presented elsewhere. In any case, that does not jeopar- 
dize our detection and identification of the quasi-molecular 
lines reported in this paper. 

The model plotted in Fig. 3 was normalized by a free pa- 
rameter in order to fit the F USE spectrum. With that nor- 
malization, our model fits as well the shape of the IUE data 
of G226-29, at longer wavelengths (~ 1250 A-2000 A). 
However, we should point out that we need to divide 
the IUE flux by a factor of <~ 1.3 to bring it to the 
model. Small uncertainties in the background subtrac- 
tion or misalignment of the target on the slit during the 
FUSE and IUE observations or both, could account for the 
flux discrepancy. It is worth mentioning that Hebrard et 
al. (2002a) fitted the FUSE and IUE data of the cool DA 
white dwarf G231— 40 without applying any correction fac- 
tors to the fluxes. Given that the flux of G231— 40 is more 
than 20 times higher than the flux of G226— 29, system- 
atic uncertainties in the background subtraction are signif- 
icantly reduced. It is also easier to quantify the alignment 
of the target on the slit as a function of time and maintain 
a reasonable photometric accuracy for a relatively bright 
target such as G231— 40. We considered that a difference 
by a factor ~ 1.3 in the flux calibration of observations 
performed with 2 different instruments of a faint target is 
acceptable. Our model with T e ff — 12000 K and log g = 7.9 
fits as well the UV spectra (~ 1250A-2500A) obtained 
with HST/FOS and presented by Kepler et al. (2000). We 
note that this HST/FOS spectrum was normalized in or- 
der to fit the IUE spectra. 

5. CONCLUSIONS 

We have presented the first far-UV spectrum of a ZZ 
Ceti star, namely G226-29, performed with FUSE. We 
fitted these data using theoretical calculations of stellar 
profiles that include perturbations by both neutral hydro- 
gen and protons. The observed spectrum is well fitted, 
allowing us to reproduce the global shape of the stellar 
continuum, together with the H~ absorption at 1130 A, 
the H^ quasi-molecular satellite 1076 A (and marginally 
the one at 1058 A), and a broad H 2 CI satellite at 1150 A. 
This is the first detection and identification of that H 2 
quasi-molecular satellite of Ly-/3. This feature has been 
predicted but never observed up to now. The quality of 
the fit allows that absorption feature to be unambiguously 
identified. 

This discovery confirms our theoretical prediction and 
gives confidence concerning the accuracy of the fundamen- 
tal molecular data we used. In combination with existing 
optical and ultraviolet data, fits using such models will be 
used in the future to allow stellar parameters (tempera- 
ture, gravity) to be constrained more accurately, and will 
provide insights crucial in the modeling of cool white dwarf 
atmospheres. 
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Fig. 1. — Line profile of Ly-/3 perturbed by neutral hydrogen and protons (solid line). The dashed line shows only the contribution from 
neutral perturbers. I w is the normalized line profile proportional to the absorption coefficient as described by Allard et al. (1999). 
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Fig. 3. — FUSE spectrum (histogram) compared with the theoretical model (solid line) with the new Ly-/3 broadening including the 
quasi-molecular satellites. The H2 quasi-molecular satellite at 1150 A is detected for the first time. The dotted line is the model with no H2 
transitions included. The emission lines (labeled ©) are due to the airglow emission in the high terrestrial atmosphere (Fcldman ct al. 2000). 



